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We present a new global QCD analysis of nuclear parton distribution functions. In addition 
to the most commonly analyzed data sets for deep inelastic scattering of charged leptons 
off nuclei and Drell Yan di-lepton production, we include also measurements for neutrino- 
nuclei scattering as well as inclusive pion production in deuteron-gold collisions. The 
emerging picture is one of consistency, where universal nuclear modification factors for 
each parton flavor reproduce the main features of all data without any significant tension 
among the different sets. 



1 Motivation 

In the last few years, significant progress has been made in obtaining nuclear PDFs (nPDFs) 
from data. In addition to the theoretical improvements routinely used in modern extractions of 
free proton PDFs, such as the consistent implementation of QCD corrections beyond the LO [I] 
and uncertainty estimates [3] , the most recent determinations of nPDFs have also extended 
the types of data sets taken into account, moving tovifards truly global QCD analyses of nuclear 
effects O m O |6] . The addition of novel hard probes to the fit does not only lead to better 
constrained sets of nPDFs and allows one to study the nuclear modification to the different 
parton species individually, but also tests the assumed process independence of nuclear effects. 

The deep inelastic scattering (DIS) of charged leptons off nuclear targets not only initiated 
all studies of nPDFs but still provides the best constraints on nuclear modifications for quark 
distributions. Upon combination with available data on Drell Yan (DY) di-lepton production 
off nuclear targets, a better discrimination between valence and sea quarks can be achieved. 
However, DIS and DY data only loosely constrain the nuclear modifications to the gluon density 
because they cover a too small range in the hard energy scale Q. To remedy this situation, 
data from BNL-RHIC for inclusive pion production in deuteron-gold (dAu) collisions have been 
included in the analysis of nPDFs performed in Ref. [3 . Not surprisingly, these data have a 
significant impact in the determination of the gluon distribution. The corresponding nuclear 
modification for gluons turned out to be much more pronounced than in previous estimates and 
also much larger than those found for all the other partonic species. 

Another promising avenue for significant improvements is neutrino induced DIS off iron and 
lead targets, with results available from NuTeV, CDHSW, and CHORUS [?]■ These data receive 
their importance from their discriminating power between nuclear modifications for quarks and 
antiquarks and have been included in a series of analyses in Refs. [H E]. Unexpectedly, the 
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correction factors obtained from neutrino scattering data are found to differ significantly from 
those extracted with charged lepton probes [4j |6]. At variance with these results, Ref. [5] 
confronts the neutrino DIS cross sections with nPDFs obtained in [3] without any refitting and 
finds no apparent disagreement. 

The novel global QCD analysis of nPDFs presented here [7j incorporates in a comprehensive 
way all of the above mentioned improvements and data sets. The resulting nPDFs at next- 
to- leading order accuracy supersede previous work presented in [I] . We adopt a contemporary 
set of free nucleon PDFs [8] as our reference distribution to quantify modifications of PDFs 
in nuclei. As in [5], we use a general mass variable flavor number scheme to treat charm and 
bottom quark contributions in our analysis. We use the Hessian method [S] to estimate the 
uncertainties of the nuclear modification factors and examine critically their range of validity. 



2 Framework 

Throughout the analysis, we make the usual assumption that theoretical expressions for mea- 
sured cross sections involving a nucleus factorize into calculable partonic hard scattering cross 
sections, identical to those used for processes involving free nucleons, and appropriate combina- 
tions of non-perturbative collinear parton densities and fragmentation functions. The nPDFs 
ff{x, Qo) at an initial scale Qo — 1 GeV are related to proton distributions ff{x, Qo) through 
a multiplicative nuclear modification factor Rf{x,Qo) as 

f^'{x,Qo) = Rfix,Qo)ff{x,Qo) , (1) 

where x is the usual DIS scaling variable for free nucleons. Both valence quark distributions 
are assigned the same nuclear modification factor R^{x, Qq) which we parametrize as 

R^ix, Ql) = £1 (1 - xf^ (1 + £2(1 - xf^-) (1 + a„(l - xf^) . (2) 

We also assume that the light sea quarks and antiquarks share the same correction factor 
Qq). No significant improvement in the quality of the fit to data is found by relaxing this 
assumption. We choose another factor R^{x,Qq) to parametrize medium effects for gluons. 
An excellent description of the data is achieved by relating both Rf and R^ to Rf specified 
in Eq. allowing only for a different normalization and modifications in the low- a; behavior. 
Hence we choose, without any loss in the quality of the fit, 

Rf{x, Ql) = Rfix, Ql) ££i±^ , Rf(^^^ Ql) = rA^:,^ qI) . (3) 

We note that the coefficients £1 and £2 in Eq. ([2]) are fixed by charge conservation, and if we 
further constrain £5 and tg to be equal, which, again, has no impact on the quality of the fit, 
£s is fixed by momentum conservation. The A dependence of the remaining free parameters 
^ € {a«, tts, Qfg, /3i, /32, /^a, Ql), fls, Og} is parametrized in the usual way [I] as ^ = 7^ -|- A^A''^ . The 
very mild A dependence found for some of the ^'s allows us to further reduce the number of 
additional parameters by setting Sa^ = Sa^ and = Sa^ , leaving a total of 25 free parameters, 
which are obtained by a standard minimization, without artificial weights for certain data 
sets, i.e. Wi = 1, and with statistical and systematic errors added in quadrature in A^: 
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3 Results 
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The total for the optimum fit was found to 
be 1544.7 for 1579 data points (x^/d.o.f. = 
0.994). All data sets are adequately re- 
produced, well within the nominal statistical 
range — n±\/2n with n the number of data. 
More specifically, the partial contribution to 

of all the charged lepton DIS data amounts 
to 897.52 units for 894 data points, for neu- 
trino DIS we find 488.20 units compared to 
532 data points, DY observables amount to 
90.72 units for 92 points, and pion production 
in dAu collisions adds another 68.26 units to 

for 61 data points. 

In Figs. 1-3 we show some examples of the 
good agreement between the fit and charged 
lepton DIS, neutrino DIS, and hadropro- 

duction data, respectively; see Ref. [7^ for Figure 1: Comparison to charged lepton DIS 
details. The remarkable agreement with 
charged lepton DIS data, shown in Fig. [1] is 




1.6 ; 

1.4 I 

1.2 \ 

I : 



- 1 

?0.8 



0.3 
0.2 
0.1 



rV] \=0.0I5- 


W 3 

r* x=0.045- 

nl 1 1 r 


HfT' x=0.080- 


] 

\_ x=U.I25J 


^'1 1 1 [ 

r fSl#4. 1 

: x=0.175 : 

nl 1 1 r 


L x=0.225 J 


r x=0.275 -. 


i'\ 1 1 [ 

- 

r x=0.350- 

"1 1 1 r 


r x=0.450t 


; A NuTeV^^*^**»»j|ii 

r □ CHORUS (Pb) 

L o CDHSW (Fe) J 


t'\ 1 1 t 

1- x=0.650-: 


1- x=0.750-: 

L this fit (Fe) J 

: this fll (Pb) 



a common feature of all nPDFs analyses. 
Neutrino DIS data for the averaged struc- 
ture function {F2^ + F^^)/2 are well repro- 
duced within the experimental uncertainties 
both in shape and in magnitude, see Fig. [51 
The only noticeable exception are the CDHSW 
data at values below 10 GeV^ where they 
exhibit a rather different slope than the other 
data. In fact, in this region it appears 
to be impossible to simultaneously fit all data 
sets equally well, suggesting some systematic 
discrepancy among the different neutrino data 
which needs to be investigated further. Data 
for the averaged structure function are also 
well described by our fit 7 . 

In general, results from dAu collisions 
are significantly less straightforward to inter- 
pret in terms of nuclear modification factors. 
Each value of px samples different fractions of 
the contributing partonic hard scattering pro- 
cesses, integrated over a range of x. Furthermore, since pr sets the magnitude for the factor- 
ization scale, the ratios reflect also the energy scale dependence of the effects. Apart from the 
nuclear modifications of parton densities, accounted for by the nPDFs, the cross sections are in 
principle also sensitive to medium induced effects in the hadronization process. 

Assuming factorizability for a given nucleus, such final-state effects can be absorbed into 
effective nuclear parton-to-hadron fragmentation functions (nFFs). The solid lines in Fig. 3 
represent the result of our best fit of nPDFs using the nFFs of Ref. [TU] . The fit follows well the 
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Figure 2: Comparison to neutrino DIS data 
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rise and fall of the ratio at small and high pT , respectively, but falls somewhat short in repro- 
ducing the enhancement found at medium pT- Owing to the large experimental uncertainties, 
the for this subset of data is nevertheless 
good, Xd^u/"' = 1.12, in particular, if com- 
pared to the outcome of an otherwise similar 
fit using vacuum FFs [TT] where X^^^/n = 
1.37. Data for 7r° yields in dAu collisions were 
first incorporated by EPS ^ and found to pro- 
vide a vital constraint on Rf^- At variance 
with our approach, the authors in [3^ disregard 
any medium modifications in the hadroniza- 
tion and assign a large weight uidAu in Eq. ^ , 
which drives their observed large nuclear mod- 
ifications of the gluon density. Our i?^" ex- 
hibits only moderate nuclear corrections. 

Uncertainties in the extraction of our 
nPDFs have been estimated with the Hes- 
sian method !£] for a tolerance criterion of 
Ax^ — 30 and found to be rather large 0, 
in particular, when compared to the present 
knowledge of free proton PDFs. As always, 
these estimates are only trustworthy in the re- 
gion constrained by data, i.e., x > 0.01. In 
particular, prompt photon and DY di-lepton 
production in dAu and pPb collisions at RHIC and the LHC, respectively, will help to further 
constrain nPDFs in the future; see, e.g., Ref. ^ for some quantitative expectations. 
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